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[T I Severe storms, tropical cyclones, and associated tornadoes, floods, lightning, and mi- 

^ crobursts threaten life and property. Reliable, precise, and accurate alerts of these phe- 

^ nomena can trigger defensive actions and preparations. However, these crucial weather 

,— I phenomena are difficult to forecast. The objective of this paper is to demonstrate the 

^ ^ potential of 4d-VAR (four-dimensional variational data assimilation) for exigent fore- 

I casting (XF) of severe storm precursors and to thereby characterize the probability of 

^ a worst-case scenario. 4d-VAR is designed to adjust the initial conditions (IC) of a nu- 

merical weather prediction model consistent with the uncertainty of the prior estimate of 
. ^ the IC while at the same time minimizing the misfit to available observations. For XF, 

^ the same approach is taken but instead of fitting observations, a measure of damage or 

4^ loss or an equivalent proxy is maximized or minimized. For example, XF of maximized 

, ^ significant tornado parameter (STP) would delineate relative probabilities of the threat of 

tomadogenesis as a function of time and place. To accomplish this will require develop- 
^ ment of a specialized cost function for 4d-VAR. When 4d-VAR solves the XF problem a 

\D by-product will be the value of the background cost function that provides a measure of 

^ the likelihood of occurrence of the forecast exigent conditions, and the value of the STP 

cost function that provides an estimate of the likelihood of tomadogenesis. 4d-VAR has 
^ been previously applied to a special case of XF in hurricane modification research. A 

2 summary of a case study of Hurricane Andrew (1992) is presented as a prototype of XF. 

^ The study of XF is expected to advance forecasting high impact weather events, re- 

fine methodologies for communicating warning and potential impacts of exigent weather 
ITpi events to a threatened population, be extensible to commercially viable products, such 

as forecasting freezes for the citrus industry, and be a useful pedagogical tool. Further, 
^ by including parameter sensitivity in the adjoint model, XF could be extended to include 

parametric uncertainty. 
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1. Introduction 

Forecasting high impact weather associated 
with severe thunderstorms and tropical cyclones 
(TCs) is a challenging intellectual problem that 
has the potential to save lives and protect prop- 
erty. High impact weather associated with se- 
vere storms and TCs includes lightning, hail, 
floods, tornadoes (Fig.[T]), and microbursts. Cur- 
rently, such small-scale weather phenomena can 
not be usefully forecast by explicit dynamic 
numerical weather prediction (NWP) methods. 
NWP solves the initial value problem for the 
fluid dynamical equations that govern the evolu- 
tion of the atmosphere. In applications of NWP 
methods, data assimilation systems account for 
the uncertainties present in the models due to 
errors in the model representation of the atmo- 
sphere and in the initial conditions (IC) due to 
errors and gaps in the observations of the atmo- 
spheric parameters. While NWP model repre- 
sentation of the fine scale structure of the at- 
mosphere is steadily improving and today's high 
resolution research models can predict realistic 
severe storms and TCs, predictions of precise 
timing and location and details of internal storm 
dynamics do not agree well with reality. For ex- 
ample, forecasts of actual tornadoes would re- 
quire ultra-high resolution that is not yet fea- 
sible, as well as more capable observation net- 
works. Instead, useful tornado forecasts focus 
on prediction of the larger-scale characteristics 
of the environment of severe storms — large ther- 
modynamic instability, vertical wind shear, etc. 
Such empirical forecasts are applicable and reli- 
able to the extent that it is the large-scale envi- 
ronment that regulates the small-scale high im- 
pact weather phenomena. 

A characteristic of models of the atmosphere 
is extreme sensitivity to small perturbations of 
the IC. This sensitivity, coupled with the uncer- 
tainty inherent in models and observations, leads 




Figure 1: The Greensburg, KS Tornado. 
Greensburg, KS was destroyed on 4 May 
2007. Photo by Van DeWald, NWS. 
From www.hprcc.unl.edu/nebraska/may4- 
2007 Greensburg-Kansas -tornado .html. 

to a loss of predictability. Consequently, it is 
possible that a small change in the IC could pro- 
duce a very different forecast. That is, a fore- 
cast of a very large specific impact on people, 
property, or nature might result from a particu- 
lar (small, plausible) IC perturbation, i.e., one 
that is consistent with the uncertainties present. 
It turns out that it is possible, using a technique 
called 4d-VAR described below, to calculate the 
IC perturbation to produce a particular result. 
For example, suppose we calculate the perturba- 
tions to maximize the potential for tomadogen- 
sis. Then, if the calculated IC perturbations are 
small compared to the IC uncertainty, we must 
not discount the possibility of tornadoes form- 
ing even if the original (i.e., unperturbed) fore- 



cast did not indicate this possibility. |Henderson 
et al.| (|2005|) call this approach "exigent" fore 



casting because of the requirement of precision 
in the forecast and urgency in the response. In 
this paper we describe the potential application 
of exigent forecasting (XF) to forecasting se- 
vere weather precursors for tornadoes (Sect. |4]). 
Because of current limitations for forecasting 
very small scales, such as severe storms and tor- 
nadoes, applications of XF must first focus on 
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larger scales, such as the strength or path of a 
TC, or on the precursors to the small-scale high 
impact phenomena. In previous work we ap- 
plied a slight variation of the XF approach to 
controlling TCs and a summary of that work is 
provided to demonstrate that XF works (Sect. [3]). 

Our implementation of XF is based on 4d- 
VAR (four-dimensional variational data assim- 
ilation). 4d-VAR is designed to adjust the IC 
of a numerical weather prediction model con- 
sistent with the uncertainty of the prior esti- 
mate of the IC while at the same time mini- 
mizing the misfit to available observations. For 
XF, the same approach is taken but instead of 
(or in addition to) fitting observations, a mea- 
sure of damage, loss, or an equivalent proxy is 
maximized or minimized (Sect. [2]). For exam- 
ple, XF of maximized significant tornado pa- 
rameter (STP, Sect. |4|a| ) should be able to pro- 
vide relative probabilities of the threat of tor- 
nadogenesis as a function of time and place. 
To accomplish this, a specialized cost function 
for 4d-VAR is needed. This cost function is 
designed to measure the likelihood and sever- 
ity of tornado outbreaks, based on the synop- 
tic and mesoscale conditions conducive to their 
formation — the cost function does not measure 
tornado damage directly. That is, we use the 
STP as a proxy that captures the key environ- 
mental controls of tomadic development. When 
4d-VAR solves the XF problem a by-product 
will be the value of the background cost func- 
tion that provides a measure of the likelihood 
of occurrence of the forecast exigent conditions, 
and the value of the STP cost function that pro- 
vides an estimate of the likelihood of tomadoge- 
nesis. XF could be a useful complement to en- 
semble forecasting because the ensemble is de- 
signed to capture the overall distribution of fu- 
ture outcomes, but not necessarily the particular 
event of consequence that is the focus of XF. 

In the future, XF is expected to have very sig- 



nificant potential for societal benefits and to be 
applicable to a wide array of problems (Sect. [5]). 
We note that operational forecasting of severe 
high-impact weather and subsequent decision- 
making is extremely time-sensitive. We com- 
ment on efficient implementation briefly in 
Sect. 4jc XF has the potential to provide spe- 
cific guidance of worst case results to focus 
the forecast and decision-making activity. The 
study of XF is expected to advance knowledge 
of forecasting high impact weather events and 
to result in refined methodologies for communi- 
cating warning and potential impacts of exigent 
weather events on a threatened population. XF 
might be extended to commercially viable prod- 
ucts, such as forecasting freezes for the citrus 
industry. In addition, XF could be extended to 
deal with parametric uncertainty (e.g. , the uncer- 
tain knowledge of the surface roughness height, 
zo). Finally, XF should prove to be a useful 
device for laboratory exercises at the university 
level. 



2. 4d-VAR methodology 

For operational weather forecasting, 4d-VAR 
finds the smallest increment at the start of 
each data assimilation period so that the per- 
turbed nonlinear solution best fits all the avail- 
able data. Mathematically, 4d-VAR blends back- 
ground (i.e., a short-term forecast) and observa- 
tions by minimizing a functional. 



(1) 



with respect to a control vector x that describes 
the atmospheric state. Here Jb and Jo mea- 
sure the misfit of the four-dimensional analysis 
(i.e., the simulation that begins with x) to the 
background and observations, respectively. 4d- 
VAR solves this complex nonlinear minimiza- 
tion problem iteratively, making use of the linear 
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adjoint of the model, linearized about the cur- 
rent nonlinear simulation. Note that is defined 
by specifying the forecast error co variance — an 
unresolved problem for the storm-scale as well 
as for TCs and other severe weather. For ex- 
ploratory experiments, available climatological 
background error covariances appropriate for 
mesoscale forecasts of similar resolution and for 
similar synoptic regimes could be used. A hy- 
brid ensemble/4d-VAR would provide improved 
covariances (Sect. |5|a| ). 

4d-VAR and XF are useful and indeed possi- 
ble because the atmosphere is chaotic and un- 
predictable. Atmospheric motions occur over 
a huge spectrum of scales ranging from the jet 
stream to the trajectory of a cloud drop. Er- 
rors at different scales are coupled by nonlin- 
ear interactions. Tiny errors inevitably present 
in the large scales very quickly result in errors 
in the position of small-scale features. At the 
same time, errors present in small scales grow 
quickly and soon become unpredictable. These 
errors then effect somewhat larger scales and so 
on. The net effect is that very small perturba- 
tions can grow quickly. 4d-VAR finds perturba- 
tions that grow in just the correct way to best fit 
the available data or to best satisfy the XF con- 
straints. 

A few changes must be made to 4d-VAR for 
XF ( [Henderson et al.|[2005| ): the control vector 
remains x, 7^ is unchanged, Jg may be present 
or not depending on whether more recent obser- 
vations are available, and 7j ("d" for "damage") 
is added to measure the "benefits" minus "costs" 
resulting from the forecast as defined by the cus- 
tomer. The functional to be minimized for XF is 
then 

J = Jt + Jo + WdJd- (2) 

The definition of should be nondimensional 
and Wd is an adjustable weight. As is in- 
creased, the value of Jj determined will de- 
crease and vice versa. Note that is defined 



to calculate a worst case scenario: minimiz- 
ing J will simultaneously maximize costs, mini- 
mize benefits, and minimize increments with re- 
spect to the background and observations. By 
contrast, in our hurricane modification experi- 
ments (Sect. [3]), by defining to measure dam- 
age alone, our solutions must minimize damage. 
That is, we calculated a "best case scenario" to 
determine perturbations to reduce damage. 

In 4d-VAR data assimilation, the final value 
of Jh + Jo objectively quantifies the likelihood 
of the calculated changes to x, based on the most 
recent background, the observations, and knowl- 
edge of instrument and model error. This is the 
case since the log of the probability of an at- 
mospheric state, given the background, the ob- 
servations, and the associated uncertainties, is 
directly related to Jh + Jo according to Bayes' 
theorem when suitable assumptions are valid 
( |Lorenc||79 86). Note that smaller values of J 
indicate higher probability. Since this applies to 
any atmospheric state, the final value of Jh in the 
XF case will still be related to the likelihood of 
the exigent solution occ uring. It is expected that 
J has a ;t^-distribution (Desroziers and Ivanov 



20011 [Muccino et al.|2004l ). In the XF situation. 



the final value of Jh will still be related to the 
likelihood of the exigent solution occuring. As 
described in the next paragraph, by calculating 
the final value of Jh as a function of the region 
used in the definition of J^, we will obtain a map 
of the relative probability of the exigent event 
(e.g., of the occurrence of tornadogenesis). 

Ideally the objective evaluation would be 
combined with a subjective evaluation based 
on forecaster experience — including pattern 
matching — that would compare the new analy- 
sis and observations to determine if the perturba- 
tion is consistent with hypothetical, but plausi- 
ble, dynamical processes or additional observa- 
tions. A series of solutions with increasing wj, 
the weight given to Jd, gives a series of solutions 
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that are increasingly unlikely {i.e., larger Jj,) and 
associated with increasing (tornado) threat. One 
would stop increasing Wd when the subjective 
evaluation indicates a nonrealistic solution or 
when the calculated likelihood becomes smaller 
than a prespecified lower value. Plotting the val- 
ues, respectively, of and Jt, at the minimiz- 
ing solution, with respect to for each location 
and forecast time, will allow us to highlight the 
locations and times of greatest threat. 



18 


















\25\ 







Figure 2: Forecast tracks of Hurricane Iniki. 
The center of the simulated hurricane is plot- 
ted for each hour for to = h to 24 h for the 
unperturbed simulation (black dots) and for 
two target experiments, one allowing pertur- 
bations to all variables (orange) and one al- 
lowing only temperature perturbations (green). 
For reference, these tracks are plotted over the 
wind speeds of the repositioned cyclone used 
as the goal. Wind speeds corresponding to 
different Saffir-Simpson categories are plotted 
with different colors: grey for tropical depres- 
sions (> 12 ms^^), green for tropical storms 
(> 17 m s^^), yellow for category 1 hurricanes 
(> 33 m s^^) and r ed for category 2 hurric anes 
(> 43 m s 1). After|Hoffman et al.|(|2006b|). 



3. An XF prototype: Weather 
control for hurricanes 

Experiments we conducted demonstrate the 
ability of 4d-VAR to calculate the influence of 
perturbation on the future path or intensity of 
a simulated TC. In "target" experiments, the 
MM5 4d-VAR system determines the optimal 
atmospheric state trajectory, which simultane- 
ously minimizes the size of the initial perturba- 
tion and the difference (using a quadratic norm) 
between the new forecast state and a goal atmo- 
spheric state in which the simulated cyclone has 
been repositioned. That is, Jd in these exper- 
iments is similar in form to Jt,, but instead of 
measuring the difference between the analysis 
and the background, Jd here measures the dif- 
ference between a short (in the case presented 
later six-hour) forecast and the goal or target. In 
experiments for Hurricane Iniki, the simulated 
hurricane successfully matched the goal state 
at the end of the 4d-VAR period, then contin- 
ued on a track parallel to its original track, and 



missed the Island of Kauai as desired ( [Hoffman 
et al.||2006bl ). Perturbations to the initial condi- 
tions are small relative to the hurricane. When 
changes are restricted to the temperature field 
alone, the technique is less successful (Fig.|2]). 

In "damage cost function" experiments, 4d- 
VAR simultaneously minimized the size of the 
initial perturbation and an estimate of property 
loss that depends on wind speed. In these exper- 
iments Jd is given by 

■^d = ~J J D(x,y,t)P(x,y), (3) 

At A 

where At is the time period of interest and A is 
the area of interest. The fractional wind damage 
D is taken to depend only on the lowest model 
layer wind speed at location {x,y) with the func- 
tional form shown in Fig. [3| and the property 
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values P{x,y) were assigned based on a land use 



data base. (Details are given by Henderson et al. 
( I2005D .) 




40 BO 
Wind Speed (m/s) 



Figure 3: Wind cost function. In our model of 
wind damage, wind speeds less than 25 ms^^ 
cause no damage and wind speeds in excess of 
90 ms^^ cause total loss. Between these two 
threshold wind speeds a cosine curve parame- 
terizes the percent of wind damage. 

In experiments for Hurricane Andrew, the 
hurricane surface winds decrease over the built- 
up area at landfall ( [Henderson et al. [2005 ). 
Additional experiments explored the ability of 
other state variables to effect desired changes 
to Hurricane Iniki. Perturbations restricted to 
winds alone, to temperature alone, or even to 
temperature only outside the center of the hurri- 



cane were found to be effective ( [Hoffman et al. 



2006a[ ). Furthermore, vertical velocity and hu- 
midity perturbations alone were ineffective at 
reducing damaging winds. An experiment us- 
ing perturbation pressure alone substantially re- 
duced the extent of damaging winds. While this 
result is similar to experiments that used solely 
temperature in the control vector, the minimiza- 
tion failed to converge to within specified nu- 
merical limits and thus should be considered 
less robust. 

The optimal perturbations usually include 



quasi-axisymmetric features centered on the 
hurricane (Fig. |4]). It appears that the pertur- 
bations then evolve as concentric wave distur- 
bances that propagate to a focus at the hurricane 
center, and convert the kinetic energy of the hur- 
ricane into thermal potential energy at the appro- 
priate time. The hurricane surface winds regen- 
erate soon thereafter, so a continuous series of 
perturbations may be needed in practice (Fig. [5]). 




Figure 4: Initial perturbation for Hurricane An- 
drew. In this color volume rendering of opti- 
mal initial temperature perturbations for Hurri- 
cane Andrew, the red and blue surfaces enclose 
the volumes where the perturbations are greater 
than 0.6 °C, and less than —0.6 °C, respectively. 
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Figure 5: Hurricane Andrew surface wind fields. Evolution of the surface wind field for the un- 
perturbed simulation and for the controlled experiment that minimized damage between 4 and 6 h 
by introducing temperatureperturbations only. Wind speed is color coded according to the Saffir- 
Simpson scale (as in Fig. [2]) and the damaging wind contour 25 m s^^ is plotted at 4, 6, and 8 h 
after the perturbation is introduced. In the experiments only wind speeds above 25 m s^ result in 
property damage. Note that at 6 h there are no damaging winds over land areas in the controlled 



experiment. After Henderson et al. (2005 1 
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4. Application to tornadogene- 
sis precursors 

A combination of inaccuracy in the initial con- 
ditions (IC), model errors, and the short pre- 
dictability time scale for severe storms all com- 
bine to make tornado forecasting very diffi- 
cult. Especially challenging are parameterizing 
the moist microphysics and estimating the IC 
for water vapor, cloud, and precipitation since 
in stormy environments these have small spa- 
tial scales (e.g., [Weiss et"aL]|2007| ). As a re- 
sult, prediction of severe weather generally uses 
an ingredients-based forecast method wherein 
coincidence of multiple individual atmospheric 
predictors increases forecaster confidence in the 
likelihood of severe weather. Historically, the 
risk of severe weather is maximized when and 
where the ingredients coincide. The discussion 
here will focus on one aspect of this overall 
problem — the IC uncertainty for dry-line tor- 
nadogenesis, where the empirical ingredients- 
based approaches work best. The XF approach 
has greater applicability than this one focus. 
In other parts of this paper, we sketch the ex- 
tensions needed to handle directly forecasting 
storm scale high impact weather elements as 
well as other high impact or high value weather 
events such as landfalling hurricanes and freez- 
ing conditions in citrus growing regions. 

Ingredients-based indices, due to their nature, 
tend to focus on prediction of sensible weather 
conditions, such as the likelihood of supercell 
thunderstorms or significant tornadoes, whereas 
individual indices or ingredients describe more 
esoteric atmospheric features, such as the mag- 
nitude of thermal instability. 

a. Cost Junction for dryline tornadogenesis 

XF of tornadogenesis based on 4d-VAR requires 
a cost function based on an ingredients-based 



index. One ingredients-based index that can 
define J^t is the significant tornado parameter 
(hereafter, STP). The STP combines individ- 
ual measures of atmospheric thermal stability, 
boundary layer moisture and horizontal wind 
fields. This index has been shown subjectively 
to discriminate between convective events in- 
volving discrete supercells that produce signifi- 
cant tornadoes and those that produce weak tor- 
nadoes or none at all. The requirement that 
the mode of the expected convection be discrete 
supercells limits the applicability of STP-based 
XF. In other situations, other more appropriate 
indices should be used. 

As defined by [Thompson et al.] ( |2003| ), 



STP 



MLCAPE 



SHR 



lOOOJkg-i 20ms-i 



SRH 



X 



2000 - MLLCL 



100 m^ s-2 
150-MLCIN 



1500 m 



X 



125 J kg 



-1 



(4) 



where MLCAPE is the mixed-layer convective 
available potential energy (CAPE), SHR is the 
magnitude of the 0-6 km vertical shear of the 
horizontal wind, SRH is the magnitude of the 
0-1 km storm relative helicity, MLLCL is the 
mixed-layer lifted condensation level (LCL), 
and MLCIN is the mixed-layer convective in- 
hibition (CIN). Values of STP greater than one 
are associated with the potential for significant 
tornadoes. Each of the components of STP has 
a history of use in tornado prediction and each 
describes characteristics of the local storm envi- 
ronment as detailed in Table [U We believe that 
use of the STP as a penalty function may be es- 
pecially relevant to characterizing the likelihood 
of CIN erosion — a frequent forecasting problem 
and one that strongly modulates the overall ex- 
tent of convection. The definition of STP is em- 



pirical and continues to be refined by Thompson 



et al. (2004 2007). However, for a preliminary 
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study the definition given in Eq. (|4]) should be 
adequate. A straightforward approach to defin- 
ing Jd so that minimizing J defined in Eq. (|2]) 
will result in maximizing the STP is to average 
STP over a target area A, and time interval At: 



1 1 



STP. 



(5) 



Af A 



Typically A will be a county of interest, and At 
a period of interest lasting from one to several 
hours. 

As defined by Eq. (|5]), Jd is a highly non- 
linear function of the atmospheric state vector, 
since not only is it a product of the five separate 
terms in Eq. (|4]), but the individual terms them- 
selves are highly nonlinear: CAPE and CIN 
correspond to the areas between a lifted parcel 
and the environmental sounding on a thermo- 
dynamic diagram, and are highly sensitive to 
changes in the mixed layer quantities that de- 
fine the parcel. Similarly, the storm-relative he- 
licity depends on an estimate of storm motion, 
which in turn depends on the depth of the con- 
vectively unstable layer. Formally, the nonlin- 
earities in can be treated in the same way 
4d-VAR handles nonlinearities in model param- 
eterizations and in complex observation opera- 
tors. In the framework of iterative 4d-VAR, the 
full nonlinear formulation is used in the outer 
loop, and linearized formulations (forward and 
adjoint) are used in the inner loop ( |Lorenc|1997 



Rabier et al.|2000| ). However, strong nonlineari- 
ties, such as "on-off" switches from stable to un- 
stable conditions, can lead to convergence prob- 
lems during the minimization. Therefore, to en- 
hance the performance of 4d-VAR, some or all 
of the individual terms of STP should be mod- 
ified. The modified versions will be better be- 
haved in terms of being differentiable and less 
nonlinear. A general approach is to smooth out 
non-differentiable points. Note that since STP 
has been defined and refined empirically, rela- 



tively small changes to the formulation of STP 
should be acceptable for current purposes. 

b. Key dryline tornadogenesis cases 

For demonstration, the 4d-VAR XF method 
should be applied to Major Tornado Outbreak 
Days ( [Schneider et al.|[2004| ), which have his- 
torically accounted for nearly 50% of all fatal- 
ities. These are calendar days during which 
there are at least six significant tornadoes, i.e.. 



F2 and higher rating on the Fujita scale ( ,Fujita 
1971[). (It should be noted that the Fujita scale 



is used here and not the newer (February 2007) 
Enhanced Fujita Scale.) Typically, Major Tor- 
nado Outbreak Days exhibit robust signatures 
for tornadogenesis in soundings and NWP fields 
and are identified by SPC forecasts as having a 
"moderate" or "high" risk for severe weather. 
Furthermore, initial case studies should be fo- 
cused on afternoon-evening springtime convec- 
tive events in the form of discrete supercells. 
These conditions often are met near the dryline 
in the Central and Southern Plains of the United 
States. Finally, because small-scale boundaries 
in the domain caused by previous convection 
(outflow boundaries, etc.) will not be well repre- 
sented by a model with currently available reso- 
lution and parameterizations, only events which 
are substantially separate, in space and time, 
from earlier convection should be considered. 
Candidate demonstration cases are listed in Ta- 
ble [H The characteristics of the 1985 outbreak 
make it the archetypical example for our pur- 
poses. 
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Table 1 : STP ingredients 



CAPE is a measure of the instability through the depth of the atmosphere, i.e., 
Ph the positive buoyant energy available, and is generally proportional to the updraft 
^ strength of a storm. Values of CAPE over 4000 kg^^ indicate extreme instability 

likely associated with high lapse rates. 



Pi 

C/3 



The magnitude of the deep-layer wind shear, along with the amount of instability, 
regulates the overall dynamical organization and persistence of a storm. Values of 
wind shear greater than 20 m s^^ in NCEP Rapid Update Cycle Version 2 model 
soundings have been shown by [Thompson et al.| ( |2003| ) to strongly support devel- 
opment of supercells. 



C/2 



The storm-relative helicity is a measure of the potential of low-level cyclonic up- 
draft rotation for right-moving supercells, i.e., supercells that move slightly to the 
right of the mean flow for a cyclonically curved hodograph. Values over 100 m^ s~^ 



are often seen with tomadic supercells. Weisman and Rotunno (2000 1 describe how 
wind shear and helicity determine supercell dynamics. 



Lower values of environmental LCL, related to high mixed-layer relative humidi- 
^ ties, tend to increase the potential buoyancy of rear flank downdrafts in storms ca- 
^ pable of supporting significant tornadoes ( |Markowski et al.|2003] ). 



CIN is a measure of the work required to lift air to its level of free convection. 
It represents the "negative" area on a sounding and is often likened to a physical 
"cap" to explosive thunderstorm development. A number of atmospheric processes 
can erode this negative area, including synoptic scale upward motion, upward mo- 
tion along a boundary and surface heating. Values of CIN associated with convec- 
tion later in the day are often greater (closer to zero) than -250 J kg^^ This term 
strongly regulates the spatial extent of significant STP. 
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Table 2: Examples of dryline tornadogenesis 



so 

ON 
CS 



OS 
OS 
ON 



o 
o 



41 tornadoes, centered in NW Pennsylvania and north of Toronto, Ontario, formed 
ahead of a cold front during the afternoon and evening hours, killing 88 people. The 



only F5 tornado in Pennsylvania history destroyed the town of Wheatland. Forbes 
( |1985 ) and Ferguson et al. ([1986) have documented the statistics of the outbreak. 



while [Farrell and Carlson] ( [1989^ attributed an elevated mixed layer (EML) as be- 
ing a major factor in the generation of the large latent instability available to the 
convection. EMLs ( [Lanicci and Wamer||1991| ) are frequently associated with con- 
vective outbreaks in the Plains of the US, yet rarely are displaced this far to the 
northeast. This case appears to be an excellent case to study due to its overall ro- 
bust and spatially distinct synoptic characteristics, predominant supercellular mode 
of convection, and the removal of the "cap" below the EML, i.e., there is no CIN 
below the EML. 

This outbreak of 66 tornadoes in Oklahoma and Kansas highlighted the risk to 
densely populated urban areas when an F5 tornado passed through Moore, OK, a 



suburb of Oklahoma City. A total of 46 people were killed by the outbreak. Roebber 



et al. (2002) and Edwards et al. (20021) describe the subtleties in the synoptic and 
mesocale environment of the outbreak which lowered confidence in many aspects 
of the forecasting of this outbreak. 

4-5 May 2007 - On back-to-back days there were major tornado outbreaks in the 



Plains, including the first EF5-rated tornado which destroyed Greensburg, KS (Mc 



Carthy et al. 2007). Synoptic conditions were conducive to supercellular storms 
during both afternoons and evenings in approximately the same geographical area. 
On the 4th, great instability was present, but weak convergence along the dryline, 
weak forcing aloft, and a capping inversion combined to limit the number of storms 
that formed, especially farther south ( Weiss et al.|2007| ). 
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c. Mode of operation and validation 

For demonstration purposes we would apply this 
procedure to a few cases selected for study. For 
the historical cases we can examine J^i for times 
and locations surrounding the actual tornadoes 
that did occur. For validation we will then sub- 
jectively compare the maps of to the tornado 
reports. We anticipate that high values of Jd 
will be correlated with the occurrence and in- 
tensities of observed tornadoes. Such a finding 
will be necessary to consider the method use- 
ful. While forecasts of severe storms are ex- 
pected to be improved with the most sophisti- 
cated physical parameterizations and the highest 
possible resolution, for forecasting precursor sit- 
uations predominantly forced by synoptic scale 
advection, the forecast model requirements are 
much less stringent. A typical 4d-VAR setup 
would use the Weather Research and Forecast 
(WRF) model configured with the YSU (Yonsei 
University) PBL parameterization CHong et al 
2006| ) and the Anthes Kuo convection scheme, a 
single-moment microphysics scheme with mul- 
tiple liquid and ice water species, and the RRTM 
radiation package. The initial and boundary 
conditions for 4d-VAR case studies are available 
eight-times daily at 32-km resolution from the 
North American Regional Reanalysis Project 
HMesinger et al.|2006[ ). 

Calculating the worst case for different ar- 
eas and times will provide a measure of threat. 
Skilled on-duty forecasters now issue outlooks 
for medium- and high-risk convective situations 
one to three days in advance. Alternatively, the 
Short Range Ensemble Forecast (SREF) system 
can be used to identify times and areas with non- 
zero probability of values of ST? exceeding a 



critical threshold ( [Bright et aL]|2008| ). In either 
case, detailed exigent forecasts would be made 
3-12 hours in advance for the areas and times 
at risk. Typically, in Eq. (|5]) the area A will be a 
county, the time interval A? will be one hour cen- 



tered on a desired forecast time, and J defined in 
Eq. ([2]) will be minimized for all areas and time 
intervals that are considered at risk or are in the 
neighborhood of areas and time intervals at risk. 
This will produce a sequence of maps of (at 
constant 7^ as described in Sect. |2]) that we an- 
ticipate will prove very useful to the operational 
forecaster. 

The methodology described here (and in 
Sect. [2]) may be quite computer intensive, yet op- 
erations in a "warn on forecast" setting are quite 
time sensitive. There are probably many possi- 
ble approaches to increase efficiency. First, 4d- 
VAR can find the minimum faster if it starts with 
a good estimate. Having found the minimum for 
one location and time, that solution may prove to 
be useful for starting the search for neighboring 
locations and times. Second, resources should 
be concentrated on locations and times most at 
risk. For example if we first minimize J for a 
large region A and long time interval A?, then 
we should be able to identify, in that minimiz- 
ing solution, smaller regions and time intervals 
most at risk {i.e., where STP is large) for further 
analysis. And this process would then be iter- 
ated down to the smallest scales required, but 
just at the most critical locations and times. 



5. Future context and outlook 

In this paper we discussed XF of the precur- 
sors of tornadoes in an attempt to avoid missed 
forecasts. We identified dry-line tornadogene- 
sis precursors as a relevant but doable high im- 
pact weather phenomena. The concept of XF is 
new and can be extended in many ways. In the 
future, with further developments of WFR 4d- 
VAR, new classes of XF experiments will be- 
come possible (Sect. |5|a[ ). Additional potential 
applications are briefly described in Sect.[5|b 
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a. Beyond precursors 

XF may benefit operational forecasters by iden- 
tifying minor changes to model initial condi- 
tions that would strongly influence the weather 
in their forecast regions. Also, this technique 
might help to identify model or analysis sys- 
tem deficiencies when applied to poorly pre- 
dicted historical cases. The sensitivity of se- 
vere storms or hurricanes to slight changes in 
environmental flow, their complex dynamics, 
their relatively compact size, and their impact 
on lives and property make them interesting 
candidates for study from the XF perspective. 
But many other situations could be usefully ex- 
amined this way. Examples include local ex- 
treme weather events, such as heavy precipita- 
tion, strong winds and extreme temperatures, as 
well as situations in which a non-linear response 
in societal or economic costs results from mod- 
est changes to seemingly innocuous values of 
the meteorological variables, such as a reduction 
in surface temperatures to slightly below freez- 
ing in a region of citrus production. XF is appli- 
cable when the consequences of a weather event, 
be they financial, related to human safety, or oth- 
erwise, would be significant whether or not the 
meteorological situation is considered "signifi- 
cant". Ultimately the technique involves evalu- 
ating the reasonableness of calculated changes 
to the existing analysis obtained by presenting 
4d-VAR with the usual model background field, 
observations, and a cost function term related to 
the specific event of consequence. 

Despite these promising opportunities, a 
number of limitations now combine to con- 
strain the experiments that could be conducted. 
The critical limitations are the initial conditions, 
models, background cost function, and impact 
of nonlinearities. Because of these limitations 
we have described a potential demonstration 
project to forecast precursor situations predom- 
inantly forced by synoptic scale advection. The 



outlook for overcoming these limitations is re- 
viewed here. 

Initial conditions: Spinning up a realistic se- 
vere storm simulation is fraught with difficulty. 
The use of high resolution radar and satellite 
data to properly initialize a mesoscale model is 
an area of ongoing research. The success of 
the technique in application to actual exigent 
weather phenomena as opposed to precursors of 
such phenomena will require greatly improved 
capability to initialize the relevant structures in 
NWP models. 

Model errors: As described above, the effect 
of model errors is not accounted for by XF. In- 
cluding additional elements in the control vec- 
tor can allow for the possibility of model er- 



ror. For example, following the idea of Tremolet 
( 2005| ), perturbations can be introduced at in- 
tervals within the 4d-VAR interval. Also see 
Sect. lSlbliil below. 

Background cost Junction: Available back- 
ground error covariances may not be appropriate 
for the particular synoptic situations of interest. 
For background error covariances, XF can make 
use of an ensemble of forecasts, such as from 
a short-range ensemble forecast (SREF) system 
(e.g., |Grimit and Mass|2002l ). 

Nonlinearities: Incremental or not, 4d-VAR 
relies on linearizing the governing equations, 
and on the adjoint of the tangent linear model. 
This results in some difficulties. First, having 
the correctly coded adjoint is not a guarantee 
of success. For example, linear instabilities can 
sometimes be present in the tangent linear and 
adjoint versions while the corresponding non- 



linear parameterization is well behaved (Mah- 



fouf 1999). That is, a formally correct adjoint 
model can be computationally unstable. Fur- 
thermore, such an instability may only be ap- 
parent in a stressing situation, e.g., a cloud re- 
solving severe storm simulation or a high resolu- 
tion hurricane simulation. Second, even at short 
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forecast times, linearization errors can be large 
( [Tremolet 2004 1. In principle the brute force ap- 
proach can be used. It now appears that suffi- 
cient computational power has become available 
to apply this approach to research questions in 
a case study context using a detailed mesoscale 
model dMartin and Xue|2004D . 



by [Riviere et"aL](|2008") and[Mu et al.1([2009[) 



b. Further applications 

i. Forecast sensitivity studies Variants of XF 
4d-VAR have been used to examine forecast sen- 
sitivity. Singular vectors are commonly used to 



identify forecast sensitivity {e.g., [Gelaro et al. 
1998p . This is normally a linear sensitivity 



in which the moist processes are often simpli- 
fied or neglected, although some authors {e.g., 
Reynolds and Rosmond[ 2003 [ ) have examined 
how singular vectors evolve nonlinearly once 
they are determined. The choice of the norm 
that measures the size of the singular vector at 
the end of the evolution period is analogous to 
our 7^. For example, singular vectors can be 



targeted to particular regions {e.g., [Puri et aL 
2001[). To do so, the amplification during the 



forecast of a small initial perturbation is maxi- 
mized with respect some norm, such as kinetic 
energy within a particular region. The singular 
vector pattern in the initial conditions may then 
be used to design a targetted observing strategy. 

XF may be considered to be a nonlinear form 
of singular vector analysis. XF might provide 
further refinement in any study or application 
that makes use of singular vectors. For example 
the XF approach could be used to diagnose how 
and why a critical forecast failed. For this pur- 
pose J^i would measure the difference between 
forecast and key elements of what in fact hap- 
pened. The control vector could be restricted 
to small regions or just a subset of the prog- 
nostic variables to establish what part of the ini- 
tial conditions were to blame for the poor fore- 
cast. Recent work along these lines is reported 



ii. Parameter sensitivity The existing MM5 
and WRF 4d-VAR systems are geared towards 
data assimilation, and the control vector used in 
the optimization process contains just the vari- 
ables defining the initial state of the model. 
It is also of interest to consider the sensitiv- 
ity to other parts of the modeling system con- 
sidered constant during the data assimilation 
process. For example, tunable parameters of 
physics packages, or "fixed" external data such 
as the roughness length {zq) or the antecedent 
soil moisture, are not known exactly, and could 
thus also be included in the control vector and 
adjusted alongside the initial conditions. Such 
approaches have been used in oceanic {e.g., 
Smedstad and 0'Brien[[T99ll [Sheinbaum[ [T^95] ) 
and atmospheric (g.g., [Wang et al.||2005|[Norris 



and da Silva[[2007[ ) data analysis and assimila- 
tion. Within XF, this would address the possi- 
bility that plausible changes in model parame- 
ters could, at least in part, lead to a model so- 
lution that maximizes tornado genesis potential. 
The associated changes to the software would 
be isolated to two separate parts of the 4d-VAR 
system: transformations from the control vector 
used in the optimization to model state variables 
and parameters, and the physical packages af- 
fected by the tunable parameters or external data 
under consideration. 

iii. Teaching tool There are many potential 
uses of the WRF system for teaching undergrad- 
uate and graduate meteorology. The XF ap- 
proach allows one to pose "what-if" lab exer- 
cises. For example, what is the minimum and 
maximum QPF to be expected from a particular 
synoptic situation given the uncertainty in IC? 
Or how much would SST have to be reduced to 
change a category 4 hurricane into a category 3 
hurricane? 
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